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Objectives 






• To enhance incompressible flow simulation capability for developing 
Jspice *S «ped«l^. u,«.e«dy f l.« pheno-enc 

associated with high speed turbo pump. 




Shuttle advanced 

Engine phase II + redesign | 


§PennState^^ 
artificial heart 


Ventri- 


Device 




LVAD 



new design 


Jr/UffgsSBainACm 


Current Challenges 


• Challenges where improvements are needed 

- Time-rntegration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• AS the computing resources changed to parallel and dijributed 
datf orL c^mpSter science aspects become important. 

- Scalability (algorithmic & implementation) 

- Portability, transparent coding, etc. 

- Hi5h-Wdi5j1ta likely to be performed using ’super 

node" which is largely based on parallel architecture 
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Time Accurate Formulation 


/ 


^^SKsStsmctCMr- 




Time-integration scheme 

Artificial Compressibility Formulation 

• Introduce a pseudo-time level and artificial compressibility 
• Iterate the equations in pseudo-time for each time step un i 
incompressibility condition is satisfied. 

Pressure Projection Method 

■ Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 




IN53D - Incompressible N-S Solver 



** Parallel version is based on INS30-UP • 

>AAPI and MLP parallel versions 

► Structured, overset grid orientation 

► Moving grid capability 

. Based on method of artificial compressibility 

» Both steady-state and time-accurate formulations 

. and S^h-order flux difference splitting for convective terms 

. Central differencing for viscous terms 

. One- and two-equations turbulence models 

. Several linear solvers : SMRES, (SS line-relcwation LU-SSS, 

GS point relaxation, ILU(0)..,.. 

•HISTORY ,, , 

** 1982-1987 Original version of INS3D - Kwak, Chong 
** 1988-1999 Three different versions were devoped : 
INS3D-UP / Rogers, Kiris, Kwak 
INS3D-LU / Yoon, Kwak 
INS3D-FS / Rosenfeld, Kiris, Kwak 
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rim» (••c) per Iteration 
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Gen Turbopump 


^/l^BtsmiCa^ir — — ^ 
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RLV 2'^‘^ Gen Turbopump 
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- 34.3 Million Points 
H 1 .4 -Three impeller rotations 
j are completed. 

d_ _ -One complete rotation 

requires less then 3.5 days 
by using 128 CPUs on SSI 
Origin 3000. When 512 CPUs 
are utilized one rotation can 
f be completed less then 1.5 
J days. In 1999, one impeller 

i rotation would take 42 days 

by using 32 CPUs on SGI 
Origin 2000 platform. 


^ M ffr- 
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Scripting Capability 




SCRIPT DEVELOPMENT FOR TURBOPUMP SMULATIONS 


Motivation 

Significant user’s effort needed 
In complex process from geometry *9* 
to flow solver 

Objective 

Develop script system to i 

- generate grids 

- create domain connectivity Input 

- create flow solver Input 

for different components anruscr 

automatically (QgQiyiCTBYy— Turbopump Scripts} 


[ DCFj— WFLOW SOLVER 


Approach 


I grl^l I domain connectfvfty Inputj [flow aoivor Input | 


Develop one script for each component with ring Interface 
between components => easy plug In for different designs 
and combinations of components 




Scripting Capability 


SCRIPT GENERATION 

Disadvantages 

> Require expertise to build scripts the first time 
Advantages 

> Allow rapid re-run of entire process 

> Easy to do grid refinement and parameter studies 

> Easy to try different gridding strategies 

> Documentation of gridding procedure 
Td scripting language 

> Works on UNIX, LINUX and WINDOWS 

> Integer and floating point arithmetic capability 

> Modular procedure calls 

> Easy to add GUI later If needed 
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Scripting Capability 


fmHimKllCftitf- 


INPUT AND OUTPUT 


> profile curves for hub and shroud In PLO'np format 
(rotated by script to form surface of revolution) 

> blade and tip surfaces In PLOT3D format 

> Parameters that can be changed 
- number of blades and sections 

_ global surface grid spacing As (on smooth r®glorw) 

? ^ nHonAfiaAnt BnO S 


- giODBI aUMBUV yiiw ■ Wnrrwk 

_ local surface grid spacing, some Independent and *on^ 
expiUsed as multlp^ of As (leading/lralling jrtpa, etc.) 

- normal wall grid spacing (viscous, wall function) 

- marching distance 
-grid stretching ratio 


Output 


> overset surface and volume grids for hub, shroud, blades 

> object X-rays for hole cutters nw n 

> domain connectivity namelist Input for OVERFLOW-D 


Scripting Capability 

^/^S 8 S 98 reH Cealer 

INLET GUIDE VANES (N repeated blades, no tip clearance) 












Manual Script (fine) Script (coarse) 

No. of pts (million) 7.1 5.8 

1 1 --|. time * 1 day 43 sec. 20 sec. 

rnSHi gISSetry deJ. to DCF Input Wllh SGI R12k 300MHz CPU) 
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Consortium Impsllsr-Dif fuser 


/A 

-Geometry and operating conditions 
(M. Williams @ Boeing , and D. Domey NASA-MSFC) 
-Impeller : 6 long blades/ 6 short bla des^ 

-Diffuser : 13 blades 





x- 

^fm^stsiAUefgt 


-R„l.do«al speed :6322 RPM / M.^ Flow 
Re • 1 37x107 / = 4.5225 inches, - 249.5 tt/sec 

-39 Million grid points / 41 zones in overset grid systems 
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STATIC/DYNAMIC STRESS ANALYSIS 
FOR TURBOPUMP SUB-SYSTEMS 


NEW CONDITIONS 



INTERFACE 


STRUCTURAL 


ZiFFER/GRID 

▲ 

► 

LOADS 


CFD GRID 
PRESSURE 
TEMPERATURE 

I 


FEM GRID 

I 


FLUIDS STRUCTURES 
INS3D NASTRAN/ANSYS 


STRESS I 
ANALYSIS! 


Summary 




•Unsteady flow simulations for RLV Sen baseline turbopump for three 
impeller rotations are completed by using 34.3 Million grid points model. 

•MPI/OpenMP hybrid parallelism and MLP shared memory parallelism has 
been implemented in INS 3D, and benchmarked. 

•For RLV turbopump simulations more than 30 times speed-up has been 
obtained. 

•Moving boundary capability is obtained by using DCF module. 

• Scripting capability from CAD geometry to solution is developed. 

•Unsteady flow simulations for advanced consortium impeller/diffuser by 
using 39 Million grid points model are currently underway. 1.2 impeller 
rotations are completed. 

•Fluid/Structure coupling is initiated. 
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